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Abstract 

In the present paper fractional Hamilton-Jacobi equation has been derived for dynamical systems in- 
V^ . volving Caputo derivative. Fractional Poisson-bracket is introduced. Further Hamilton's canonical equa- 

^Nj ' tions are formulated and quantum wave equation corresponds to the fractional Hamilton-Jacobi equation 

is suggested. Illustrative examples have been worked out to explain the formalism. 
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c^ . 1 Introduction 

In seminal papers Riewe [1,2] has formulated Lagrangian and Hamiltonian mechanics to include derivatives of frac- 
tional order [3-9]. It has been shown that Lagrangian involving fractional time derivatives leads to equations of 
motion with non conservative classical forces such as friction [1]. Motivated by this approach many researchers have 
^ \ explored this area giving new insight into this problem [10-16]. Agrawal [12] has developed fractional calculus of 

variations dealing with problems in which either the objective functional or the constraint equations or both contain 
at least one fractional derivative term. Agrawal [12] has dealt with Lagrangian involving Riemann-Liouville (R-L) 
fractional derivatives. R-L derivatives are nonlocal. R-L derivative of a constant is not zero, and in many applica- 
tions it involves fractional initial conditions which are nonphysical. For these reasons Caputo derivative [8, 9] has 
r*~" ' widely been used in recent literature. Agrawal [13, 14] in the recent papers has presented fractional Euler-Lagrange 

equations involving Caputo derivatives. In conclusion it is emphasized that both (the R-L and Caputo) fractional 
derivatives arise in the formulation, even when the fractional variational problem is defined only in terms of one type 
of derivative. Thus fractional boundary conditions may be necessary even when the problem is defined in terms of 
Caputo derivative. Further fractional Hamiltonian formulation has been developed in terms of Caputo derivatives 
L^ " by Baleanu and coworkers [15,16]. As a pursuit of this in the present paper we investigate the fractional Hamilto- 

rS , nian involving Caputo derivative and derive the Hamilton-Jacobi equations. Poisson brackets constitute important 

jrt ' part of Hamiltonian mechanics. Entire Hamiltonian mechanics can be restated in terms of Poisson-bracket. In view 

of this a generalization of Poisson-bracket (fractional version) is suggested. Hamilton's canonical equations (frac- 
tional case) have been expressed in terms of fractional Poisson bracket. Further fractional quantum wave equation is 
suggested. Illustrative examples are presented to explain the formalism. 

2 Fractional Calculus 

Fractional calculus deals with generalizations of integer order derivatives integrals to arbitrary order. In this section 
we present basic definitions and properties which will be used in the subsequent sections [3-7]. 

Definition 2.1. // f{x) e C[a,b] and a > then [4-6] 

are called as the left sided and the right sided Riemann-Liouville fractional integral of order a, respectively. 
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Definition 2.2. Let n — 1 < a < n, then 

"^^"^(^) - rX^ {iy[ (.-/)-l^.^ ^^' (2.3) 

are called as the left sided and the right sided Riemann-Liouvtlle fractional derivative of order a respectively whenever 
the RHS exists. 

Definition 2.3. Let f{x) e C"[a, b] and n - 1 < a < n, then [4, 6] 

^D^fix) ^air"D"fix) = ^^J^^^ Ijx^tr—' (^^y f(t)dt, a<x<b (2.5) 

'^DUix) =. ir"i-Drf{x) = ^^J_^^ f\t - xr—' (^-^) " /wdt, a<x<b (2.6) 

are called as the left sided and the right sided Caputo fractional derivatives of order a respectively whenever the RHS 
exists. 



Properties: [14] 
i) 

ii) 
iii) 

iv) 

V) 



aDnf{t)+9{t))^ '^Dmt)+ ^A"ff(i). (2.7) 

aDtC = 0, c is constant. (2.8) 



/ [^D?f{t)]git)dt^ f f{t)[?D^g{t)]dt. (2.9) 

J a J a 

^DUt'af = ^^^^^^^{t-af-" (/3>a). (2.10) 

./r „A"x(t) = x{t) - J2 f(^Vr^(' ~ ''^"~'- (2-^2^ 

tl^ tD^xit) = x(t) - J2 r(a +'w) (^ ~ '^"''- (2-^^^ 



n-1 



{D^x)(a) 



J: ^D?xit) = x{t) - J2 r( +1) ^' - ")'• (2.14) 



3 = 
n-1 



,/," f i3?:c(t) = x{t) J2 %^1#^(^ - ^)'- (2.15) 

3 Fractional Mechanics 

Agrawal and coworkers [13, 16] have presented Euler-Lagrange equations for fractional variational problems defined 
in terms of R-L and Caputo derivatives. In the following section we state a theorem regarding Lagrangian involving 
left and right Caputo derivatives, which will be used in further discussion. 

Theorem 3.1. Let J[q] be afunctional of the form 

J[q] = f L{t,q, ''^D'^q, ?D^,q)dt, (3.1) 

where < a, /3 < 1 and is defined on the set of functions fix) which have continuous left Caputo fractional derivative 
(LCFD) of order a and right Caputo fractional derivative (RCFD) of order (5 in [a,b]. A necessary condition for 
J[q\ to have an extremum for a given function q{t) is that q{t) satisfies the generalized Euler-Lagrange equation: 



and the transversality conditions: 



r,a-i dL n/3-i dL 



n{t)t^Q, (3.3) 



where tD"~ denotes the fractional integral of order 1 — a. See [13] for a proof. 

4 Fractional Canonical Transformations and Generating Functions 

In this section we present the Hamiltonian formulation involving Caputo fractional derivatives. Consider the fractional 
Lagrangian given in equation (|3.ip . Then the canonical momenta pa and pp are 

dL dL .^^. 

where Pa and pp are independent. The fractional canonical Hamiltonian is 

H^p^ ^D?q + pp ?D^q - L. (4.2) 

Taking total differential of (14. 2|) and using (|4.H) . we obtain 



"^ a?"" dt 



dH = dp^ '^D^q + dp0 'rD^q - —dq - — dt. (4.3) 

Taking into account the fractional Euler-Lagrange equations (|3.2|l we get 



dH = dp, ^D^q + dpp ?D^q + ( tDb>„ + aDfpf,)dq - — dt. (4.4) 



aDtq + dpi3 ?D'^q + ( tD°p„ + aD'^Pi3)dq - — 

Equation (14. 4|) shows that i? is a function of pa, p^g, g and t. Comparing total differential of H equation (|4.4|l we 
have: [16] 

dH dL dH c r~,a dH c t^i3 dH p ^i ,. r\ 

^ = -W d^= -''^'' Wr''"' -^= ^D,p,+ .D,p^. (4.5) 

Transformation of q, Pa, p/3 into new variables Q{q,Pa,Pi3,t), Pa{q,Pa,Pp,t), F'i3{<l,Pa,Pi3,t) is canonical if there 
exists a new Hamiltonian 7i(Q,Pa,Pf),t) which satisfies modified Hamilton principle: 

5 f {Pc ^DTQ + Pp ?D^Q ~ H)dt = 0. (4.6) 

J a 

As g, Pa, and pfj are canonically conjugate, we have 

5 / (Pa :r A"g + P/3 f of g - n)dt = 0. (4.7) 

For these equations to hold, the integrands must differ by a total time derivative of an arbitrary function G, hence 

(P, f A"Q + Pp ?D^Q - H)dt - (pa ^Dl'q + pf3 ?D^q - H)dt = dG. (4.8) 

Since G is not varied at the end points, we get 

rb 



S f ^dt = S[G{b) - G{a)] = 0. 

Ja dt 



(4.9) 



The function G, which completely determines the transformation is called as a generating function. For mechanics 
involving fractional derivatives, we introduce variables qa, qp, Qa, Q/3 satisfying 

^= "^*'^' ^^ '^""^ ^r^ '"^'^' ^= '^o^- (^-^o) 

For integer order derivatives, these new coordinates are the same as the usual canonical coordinates. However, while 
dealing with fractional derivatives, the coordinates qa, qp, Qa, Qp will not be canonical, so all canonical expressions 
must be written in terms of the original coordinates aD'tq, tL>^q, aDtQ, tL>^Q. 



5 Canonical Transformation of The First Kind 

For a generating function G{qa,qp, Qa, Qi3,t) the transformation is 

(p. ^D?q + p0 ^Dlq ~ H)dt-{P^ ^D^Q + Pp ?D^Q ~n)dt = dG{q^,qp,Q^,Qp,t). (5.1) 

We have 

,^ dG , dG , dG ,^ dG , == dG , 

dG = -TT^dq^ + -rr^dqn + Tr^dQc + T^T^dQp + —-dt. 5.2 

oqa oqp dQa dQp at 

Eqs. (|4.9l5.ip yield 

dG dG dG dG p dG ^ „ ,,„, 

dqa dqp dQa dQi3 dt 

6 Canonical Transformation of The Second Kind 

Let <S be a generating function dependent on q^, qp, Pa, Pp, t. From equations (|4.10|l . (jS.lf) we have 

dG = Padqa — PadQa + ppdqp — PjsdQp + {H — H)dt (6.1) 

= Padqa - d{PaQa) + QadPa+Pl3dq0 - d{Pf)Qp) + QpdPp + {H- H)dt. (6.2) 

It is easy to observe that 

d(G + PaQa+PtiQl3)= Padqa + Q adPa + Pfidq^i + Q (idPfi + {H - H)dt. (6.3) 

Let 5 = G + PaQa + PpQii then 

d5 = Pa dqa + QadPc + ppdqp + Q/sdP^ + (H - H)dt. (6.4) 

Since 5 is a function oi qa, qi3. Pa, -P/3, i we can write 



dS _ dS _ dS _^ dS _^ dS 

r,_ — Pa, „_ — Pp, „„ — Qa, „„ — y^g, „ 

dqa dqp dPa dPp dt 



'a, Tn-=P/3. 7TB- = <5"' 7r^ = Q'3' ijr=^--^- (6-5) 



7 Fractional Poisson Bracket 

Hamiltonian mechanics can be written in terms of Poisson brackets. In the present section a generalization of Poisson 
bracket has been introduced, which is useful for generalizing fractional mechanics involving Caputo derivatives. 

Definition 7.1. If the functions F[t,q,pa,Pi3) and G{t,q,Pa,Pi3) depend on the position coordinate, fractional mo- 
menta and time, fractional Poisson (FP) bracket of F and G, denoted as [F,G]fp is defined to be: 

rr.^1 dF.dG dG, dG.dF dF. ,^,, 

F,G>P = — -— + -—-— -— + -—. 7.1 

dq dpa dpp dq dpa dpp 

The following properties can be observed: 
(a) 

[F,G]fp = -[G,F]fp, 

(b) 

[Fx+F2,G]fP = [Fi,G]fp + [F2,G]fp, 

(c) 

[Fi, [F2,F-i\Fp]FP + [F2, [F-i, Fi\fp]fp + [F-A, [Fi,F2\fp]fp = Jacohi's identity, 

(d) 

rr. 1 ,dF dF , 

lF,q]FP = -{j^ + j—), 

dpa dp0 

(f) 

dF 

[F,Pa]FP = [F,Pi3]fp = — , 

(g) 

[q,q]FP = [Pa,Pa]FP = [Pa,P,a]FP = 0, 

lPa,q]FP = [Pl3,q]FP = -1. 



8 Hamilton's Canonical Equation In terms of Poisson Bracket 

Hamilton's canonical equations in terms of Poisson brackets can be expressed as follows 

[p^,H]fp = [pp,H]pp = ^{aD^pp +t D^pc,) = ~^. (8.2) 

9 Fractional Quantum Wave Equation 

As in conventional mechanics, the Hamilton- Jacobi (H-J) equation results from a canonical transformation for which 
the new variables are constant. For integer-order derivatives, such a transformation will follow automatically if the 
new Hamiltonian Ti is identically zero, since from the equations of motion we then have 

Q-^-OP-~^ (91) 

For fractional derivatives, we can derive a similar relationship by putting 

dS 

— + H{q,p^,P0,t)=O. (9.2) 

In view of (|6.5p , (|9.2p yields fractional version of H-J equation i.e 

dS ,,, dS 9<S , „ ,„ „, 

— +H{q,-—,-—,t) = 0. 9.3 

at dqa oqi3 

Hence the quantum wave equation corresponding to Hamilton- Jacobi involving fractional Caputo derivative is sug- 
gested to be [1] 

[H{q,-ih^,-ih^,t)]^ = ih^, (9.4) 

aqa oqp at 

where t/; is a wave function. 

10 Examples 

Example 1. The total energy of the fractional oscillator is given as [17] 



H = -kx^ + -m,c{ aDtxf, 0<a<l, a < t < b, (10.1) 



The generalized energy of fractional oscillator in uniform electric field E will be 



Hfo = qEx + hix^ + ima( aDtxf, < a < 1, a < t < b, (10.2) 



where q the charge of oscillator, nia = 7 m and dimension of 7 is T" . In view of Eq. (|4.2p . Lagrangian for FO is 

1 1 

Lfo = yTriai aDt^)^ - -^kx^ - qEx. (10.3) 

The generalized Euler-Lagrange equation Eq. (|3.2p and the transversality condition Eq. (|3.3p yield: 

(a) -qE-kx + m^ tD^ ^D?x = 0, (b) tD^^ '^D^x\t=t = ei, (10.4) 

respectively. Eq. (|10.4p can be solved to get explicit expression for x(t), as follows: 
Let oj^ — -^ and - — = 7. The Euler-Lagrange equation then takes the form 

tDb aDtX = Lulx + 7. (10.5) 

To find the solution of Eq. pO.Sp with initial condition a;(0) = eo, we apply t/" to both sides of Eq. (110. 5p . Further 
using (f2T3)l and (fTOill fb'). we get 

^D?x{t) = Lul J^x{t) + \ {b - tr + ^{b - tr-\ (10.6) 

i (a + IJ 1 (a) 

Similarly, applying „/" to both sides of Eq. (|10.6p and using Eq. (|2.14p we obtain 

x{t)=ulal^ J^x+ „J[«[-£i-(6_t)"-i + --^_(6-i)"]-feo. (10.7) 

1 (a) i (a -I- 1) 



as (2;(0) — eo). Equation (|10.7|) can be thought of as Volterra-type integral that has composite integral operators and 
can be solved to obtain 

x{t) = {l-c.l al? tin-'i al-lT^Ab ' tr^' + ^, l.A b - tr] + eo) (10.8) 

i (a) 1 (a + 1) 

oo 

= J2i<^i JT ^in'i am-^Ab - tr-' + ^^^(& - tn + eo). (10.9) 

The generalized momentum for FO is 

9LfO C T^a /in 1n^ 

P"^ aCD'^ =mcc aDtX, (10.10) 



and generalized Hamiltonian in terms of Pa is 



1 2 ,1, 2 



Hfo = ^ pi + -kx' + qEx. (10.11) 



Hence the Hamilton's equations (Eq. (14. 5p ) become 

dHpO Pa Cr,a„ dHpO 1 , TT, ncy /1n10^ 

-^ = = aDtq, —^ = kx + qE= tDt Pa, (10.12) 

opa ma oq 

which are equivalent to generalized Euler-Lagrange equations given in Eq. (|10.4|l (a). 

Further we present Hamilton- Jacobi equation for this system. Eq. (|9.3p yields the Hamilton- Jacobi Equation. 

dSpO , 1 ,dSF0.2 , 1 , 2 , ni n /inio> 

^r + 2^^^i7^ + 2^" + '^^^ = °' (^°-^3) 

where ^ = aD"x{t). Assume that a solution to Eq. p0.13|l is the form of Sfo = Si(t) + S2{xa)- Then 

' (^f + ^fc^' + 9^^--^=A (10.14) 



where /? is a constant. Therefore 



So that 



^ = J2m„{p ~ ifca;2 - qEx), Si = -/3i. (10.15) 



Sfo = Xc\l2ma{(3 kx^ - qEx) - I3t. (10.16) 

Further we identify (3 with the new momentum coordinate Pa- Hence the new position variable Qa will be 

Qa = ^1^ = -^[x^^2m^{p -hx^~ qEx) - I3t] (10.17) 

m _ 

-t = 7, (10.18) 



l2mc(H ~ \kx'^ ~qEx) 
where 7 is a constant. Solving Eq. (|10.18|l for Xa we get 



maXa = (7 + t)\/2m„(/3- -kx'^ - qEx) (10.19) 



Further differentiating (|10.19p with respect to t, we get 



dx 1 

m^-^ = W2m„(/3--fcx2-gSx) =pc. (10.20) 



In view of Eq. (|9.3|l . we can identify /? with the Hamiltonian _ff. Then in view of (llO.llf) rn^^ — w2m(/3 — ^kx^ — qEx) 
Pa- But 

niaaDtX^P^. (10.21) 

Applying tD^ to both sides of Eq. (|10.2ip we get 

ma tD^ aD^x = tD^p^. (10.22) 

Using Eq. (|10.12p we obtain 

mc tD" ^D"x = kx + qE, (10.23) 



which is the same equation that we derive from Euler-Lagrange and Hamilton equations. 
This result can also be obtained using fractional Poisson brackets. 

[x,Hfo]fp = [x,-kx^ + qEx+ -■ma{'^Dtx)'^]Fp. 
Since the generalized momentum for FO is 

we can write, 

1 ^ 

\x, Hfo]fp = [x, -fcx" + qEx + -^]fp = ^= ^D^x, (10.25) 

2 Znia nia 

and 

1 ^ 

[Pc.,Hfo]fp = [pc, ^kx^ + g£a; + ^]fp = -kx - qE = - tD^Po,. (10.26) 

These equations are the same as the Hamilton's equations Eq. (|10.12p . Therefore, both the methods yield the following 
equation, for fractional oscillator in uniform electric field 

-kx-qE + matD'i aDtX^O. (10.27) 

Eq. (|9.4|l yields the following wave equation corresponding to fractional oscillator in the uniform electric field E 

(t^(-*^7F-)' + ^^^' + ''■^^)V' = i^^- (10-28) 

Ima oxa 2 at 

The wave function for the quantized FO is of the form 

iP[x,x^,t) = A[x,Xc.,t) ei'Sfo(--^-'), (10.29) 

where A{x,Xa,t) and SFo{x,Xa,t) are amplitude and phase respectively. Substituting Eq. p0.29p in Eq. (|10.28p we 

Separating this expression into real and imaginary parts, we get the two equations: 

1 ,dSF0.2 OSfo 1, 2 r^ H^ , dA .2 ^_, dA 



2ma axa at 2 2ma axa at 

Eq. (|10.3ip reduces to fractional Hamilton- Jacobi, when ft = 0. 

Example 2. Consider FO with dissipative force: 

F{x) = -7 f of a; (10.32) 

where dimension of 7 is MLT~^^'^^' . Using potential corresponding to this force we have Lagrangian as [1]: 



-ma[aDtx) ~ -kx ~i j-( '^2^^ 

2 ^ 2(-l): 



L = -m,( ''aDU) - l^kx^ - i ^( t D^ xf . (10.33) 



In view of Eq. (|3.2|l we obtain generalized Euler-Lagrange equation as 



The generalized momenta are 



-fca; + m^ tD^ "^D^x - i—^ aOf ? d'^ x = 0. (10.34) 

(-1)2 



Pa = ma aDtX and pjs = -i 5- ^ D? x. (10.35) 

(-l)f 



The corresponding Hamiltonian is 



H = \maCD'^xf + \kx^ -f i ^ (f d| xf. (10.36) 

2 ^ 2(-l)t 



Using Eqs. p0.35|) the Hamiltonian for this system takes the form 



1. 



P\ 



H = -^ + ^kx^ + ^-^ . (10.37) 



2ma 2 



2J7(-1) 



Using Hamilton's canonical Poisson bracket we have 

2 , -# 



[p^,H\ = [p„ ^ + -kx' + ^-^] = -kx = b^,H] = -{aD.^p, +t DtPa). (10.38) 

inia 2 2J7(— 1)^ ^ ^ 

Substituting Eqs. (|10.35p in Eq. (|10.38p we arrive at Eq. (110. 34p . 



11 Results and Conclusions 

Fractional mechanics describes both conservative and non-conservative systems. With this motivation a generalization 
of Poisson bracket is introduced. Further generalized Hamilton's canonical equations have been derived in the case 
of Lagrangian involving Caputo derivatives. Subsequently Hamilton- Jacobi equation has been derived and fractional 
quantum wave is suggested. 
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